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Abstract 

Consider a random matrix whose variance profile is random. This random matrix is ergodic in one 
channel realization if, for each column and row, the empirical distribution of the squared magnitudes 
of elements therein converges to a nonrandom distribution. In this paper, noncrossing partition theory 
is employed to derive expressions for several asymptotic eigenvalue moments (AEM) related quantities 
of a large Wishart-type random matrix HH^ when H has a random variance profile and is nonergodic 
in one channel realization. It is known the empirical eigenvalue moments of HH^ are dependent (or 
independent) on realizations of the variance profile of H when H is nonergodic (or ergodic) in one 
channel realization. For nonergodic H, the AEM can be obtained by i) deriving the expression of AEM 
in terms of the variance profile of H, and then ii) averaging the derived quantity over the ensemble of 
variance profiles. Since the AEM are independent of the variance profile if H is ergodic, the expression 
obtained in i) can also serve as the AEM formula for ergodic H when any realization of variance profile 
is available. 

Index Terms 

Random matrix, Wishart matrix, variance profile, asymptotic eigenvalue moments (AEM), noncross- 
ing partition. 
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I. Introduction 

Consider the linear vector memoryless model y = Hx + w, where x, y and w are the input vector, 
output vector and additive white Gaussian noise (AWGN), respectively, and H denotes the random channel 
matrix independent of w. Entries of the matrix H depend on the actual application, and the linear model 
is characterized by the joint distribution of entries in H. It is known that, if the elements of a sized N xK 
random matrix H are independent and identically distributed (i.i.d.) zero-mean random variables having 
a common variance l/N, then the empirical distribution of the eigenvalues of HH* converges almost 
surely (a.s.) to the Marcenko-Pastur law [1] when K, N — > oo with a finite ratio K/N . However, in many 
applications, entries of H have unequal variances, i.e. independent but non-identically distributed (i.n.d.). 
Examples include direct sequence-code division multiple access (DS-CDMA) and multicarrier (MC)- 
CDMA systems with frequency-flat/selective fading, multiaccess system with antenna diversity, and so 
forth. Moreover, in a multiple-input-multiple-output (MIMO) system with spatially correlated fading, the 
channel matrix with correlated elements can be transformed to another matrix having the same asymptotic 
eigenvalue distribution (AED) but with independent and unequal-variances components [2]-[4]. In these 
cases, explicit expressions for the AED of HH^ rarely exist. Fortunately, in many occasions, the AED 
can be characterized in a way that several performance metrics such as capacity, minimum mean-square- 
error (MMSE), multiuser efficiency, optimum output signal-to-interference-plus-noise ratio (SINR), etc., 
can be evaluated through solving nonlinear equations numerically. A number of works follow this line 
include [3], [5]-[ll], and a more complete list can be found in [12]. They all start from the celebrated 
Girko's theorem [5, Corollary 10.1.2], where the Stieltjes transform for the AED of HH* is obtained as 
an integral of the solution to a fixed-point equation depending on the variances of entries in H. 

In the literature, the asymptotic eigenvalue moments (AEM) of a random matrix are much less explored 
than AED. This is mainly because a moment based method reveals little underlying behavior and the 
combinatorial arguments involved are frequently horrible. Nevertheless, it appears unlikely to dispense 
with the method since, for example, there has been no alternative way of proving the behavior of the 
extreme eigenvalues [13]. Some application cases of AEM are exemplified below. 

• Let the moments = E{A fc }, k = 1, 2, ■ ■ ■ be available for K, N — > oo and finite K/N, where 
A is the random variable governing the eigenvalues of HH^. These moments are called the AEM 
of HH^. When the Carleman's criterion Y^k=\ m 2k^ 2k ^ = 00 holds, the moment sequence {nik} 
uniquely determines a distribution (AED of HH^) [14]; thus, one can use {m^} to characterize the 
system. Suppose that the goal is to evaluate the expectation of a certain function g of the random 
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variable A, e.g., g(X) = log 2 (l + 7A) for the ergodic capacity and g(X) = 1/(1 + 7A) for MMSE 
under the input signal-to-noise ratio (SNR) 7. The Gauss quadrature rule method [15] can be adopted 
to achieve the goal, where the expectation E{g(X)} is expressed as a linear combination of samples 
of g(X), and {rrik} is used to determine the coefficients in the combination and the points to be 
sampled. 

• The AEM is relevant to the design and analysis of a reduced-rank MMSE receiver [10], [16]— [21], 
which is able to reach the performance of a linear MMSE receiver with a lower computational cost. 

• The AEM plays an important role in free probability theory [22], which has been lately used as a 
powerful tool to analyze complicated wireless communication models [23]-[30]. 

Let Hi j denote the (i, j)-th entry of matrix H. The variance profile of a random matrix H is a matrix 
whose (i, j)-th entry is the variance of Hij. Consider the case that the variance profile of H is random as 
well. Then H is said to be ergodic in one channel realization if, for each column and row, the empirical 
distribution of \Hij\ 2 therein converges to a nonrandom distribution [11, Defintion 2.3]. It is known, 
if H is ergodic in one channel realization, the empirical eigenvalue distribution of HH* converges a.s. 
to a nonrandom limit whose Stieltjes transform can be given by [5, Corollary 10.1.2], where it can be 
seen the AED as well as AEM of HH^ are independent of any specific realization of variance profile of 
H. However, in many practical applications, H is nonergodic in one channel realization. In this paper, 
we assume that, although H is nonergodic in one channel realization, it is ergodic in the time domain. 
Specifically, we suppose that the variance profile of H is controlled by some random variables (denoted by 
0); conditioned on a certain variance profile, the conditional AEM converges asymptotically to quantities 
that are functions of realizations of 0. Since H is ergodic in the time domain, further averaging over the 
ensemble of the random vector 6 yields the unconditional AEM. 1 

Most previous works employing asymptotic results of eigenvalues of HH^ either have the results 
conditioned on a certain realization of H's variance profile (deterministic variance profile) or consider 
random variance profile but assume H is ergodic in one channel realization. Examples of nonergodic 
H are given below. A one-shot asynchronous CDMA system is considered in [8], where the empirical 
distribution of \Hij\ 2 in a column is a function of the random relative delay and received power of the 



'in some occasions, instead of performing the averaging process on AEM, it is more reasonable that we compute the average 
of the quantity of interest, e.g. spectral efficiency, coefficients of a reduced-rank MMSE receiver, etc. However, these quantities 
are generally related to 6 in such a complicated manner that it is intractable to evaluate the average. Thus, performing expectation 
over the AEM and then using it to obtain the interested quantity appears to be a feasible alternative. 
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user corresponding to that column. To find the optimal output SINR, the statistics of the random variables 
controlling the variance profile matrix are incorporated into the fixed-point equation governing the Stieltjes 
transform. In [11], a number of performance measures of an MC-CDMA system are analyzed, where a 
situation is taken into account that the frequency-selective fading channel is nonergodic in the frequency 
domain but ergodic in the time domain 2 . The performance measures conditioned on a specific variance 
profile matrix are averaged over the ensemble of the fading process to yield the ergodic performance 
measures. In [2], [4], mappings between the scattering radio environment of MIMO and a spare virtual 
channel matrix (or called coupling matrix) are established. The spare channel matrix contains many 
zero elements, where the number of nonzero entries indicates the channel degree of freedom, and the 
distribution of non-zero elements in the grid corresponds to the radio environments. It is desirable to 
construct a statistic model for the virtual channel matrix to govern the ensemble radio channels over time 
and space. 

Some related works about AEM computation with i.n.d. entries of H are summarized, and short 
comments are made to compare them with our work. Expressions for AEM of HH^ are derived in [10] 
for the cases that entries of H are zero-mean and i.n.d. and H is ergodic in one channel realization. The 
expressions therein are recursion based. They are obtained by expanding the fixed-point equation of the 
Stieltjes transform of AED as a power series. As a Stieltjes transform can be formulated as power series 
with moments as the coefficients, AEM are identified by equating the two power series expressions (one 
by expanding the fixed-point equation of the Stieltjes transform, and the other by a power series with 
moments as the coefficients) and comparing the coefficients of every degree at either side of the equality. 
In the current paper, for a particular realization of H's variance profile, we derive conditional AEM 
expressions in both direct and recursive forms. For the latter, our expression exhibits a simpler form than 
that in [10]. Recursive forms have the advantage of lower computation complexity when the moment 
order is high; however, it is inconvenient to perform averaging process on them to obtain unconditional 
AEM formulas. In [6], a fundamental observation is made that limits of Gaussian band matrices are 
operator-valued semicircular elements; thus operator-valued free probability can be used for determining 
their eigenvalue distribution. The works of [31], [32] employ the result of [6] and take an operator- 



2 In this case, the randomness in the spectrum of the channel matrix due to the realization of the fading process does not 
vanish asymptotically. Thus, the performance measures converge asymptotically to quantities that are functions of the fading 
realization, and further averaging over the fading ensemble yields the ergodic asymptotic performance measures when the fading 
is temporally ergodic. 
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valued free probability approach to calculate the limits of the eigenvalue distributions of a number of 
block matrices, including Wishart-type block matrices. In the intermediate process, the operator-valued 
moments in a recurrence relation can be obtained. By using a trick of putting H and in off-diagonal 
blocks of a 2 x 2 block matrix X (diagonal blocks are zero matrices), a Gaussian matrix X is obtained 
whose AEM (conditional) has a simple mapping with the AEM of HHt However, the AEM obtained 
in this way is also in recurrence, and it is inconvenient to extend to unconditional AEM. 

The rest of this paper is organized as follows. Section II lays down some necessary definitions. In 
Section III, expressions of AEM in both direct and recursive forms are derived given a certain realization 
of the variance profile of H. In Section IV, AEM formulas are obtained by means of the direct-form 
result in Section III, and two statistic models of nonergodic H are investigated. Applications are given 
in Section V, and this paper is concluded in Section VI. 

II. Preliminaries 

Definition 1 (Empirical Distribution): The empirical distribution of the vector [v\,V2, - ■ ■ , vn] is given 

as 

1 N 

where u(-) is the unit step function. ■ 
Consider m N x K matrix H = M S, where denotes the element- wise matrix product, and M 
and S are independent N x K random matrices. The entries of S are arbitrarily distributed zero-mean 
i.i.d. complex random variables with variance 1/N. We denote the (i, j)-th entries of H and M by iTjj 
and Mi j, respectively. 

Definition 2 (Ergodicity in One Channel Realization [11]): Let |_-J denote the closest smaller integer. 
For a given x £ [0, 1), let the empirical distribution of 

[\ M [xN\,l\ 2 -, \ M [xN\,2\ 2 -, ••• i \ M \_xN\,k\ 2 \ 

converge to F x (-) when K, N — > oo and K/N — > [3 a finite constant; for a given y € [0,1), let the 
empirical distribution of 

[\ M l,[yK\\ 2 , W 2t [yK\\ 2 , ■■■ , I^V^jl 2 ] 

converge to F y (-) as K, N — > oo and K/N — > [3. If the asymptotic empirical distributions F x and F y , 
x,y G [0,1), are a.s. nonrandom limits, then the random matrix H is called ergodic in one channel 
realization. ■ 
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Conditioned on a realization of M, the variance of Hij is equal to \Mij\ 2 /N. The variance profile of 
H given a certain realization of M is defined below. 

Definition 3 (Conditional Variance Profile): Consider a realization of random matrix M. Let 

V(i,j) := \M id \ 2 . 

For each N and K with their ratio K/N = (3, let : [0, 1) x [0, 1) -> R+ be a function given by 

vW(x,y) = V(i,j), l -jl< x< ^, Lzl< y< j-. (1) 

Assume that the sequence {v( N \x,y)}J$ =1 converges to a limiting bounded function v(x,y). Then, 
v(x,y) is referred to as the conditional variance profile of H given M. Furthermore, for a given 1 < 
k < K, let be a function defined in [0, 1) such that 



vF>(z) = V{i,k), ^ r <x<-. (2) 

The limit of denoted as v^[x), is called the conditional variance profile of the A;-th column of H 
given M. ■ 
By the Girko's theorem [5, Corollary 10.1.2], conditioned on M, the empirical eigenvalue distribution 
of Htf converges a.s. to a nonrandom limit whose Stieltjes transform S[z) can be expressed as 

S(z) = / u(x,z)dx, (3) 
Jo 

where u(x, z) satisfies the equation 

u(x,z) = : ——. — r . (4) 

^ M JO l+zf L u(w,z)v(w,y)dw y 

It is readily seen from (3) and (4) that, when H is ergodic in one channel realization, the conditional 
AED (as well as the quantities of our interest, AEM) of HH^ given M is invariant to the conditional 
variance profile v(x,y) and hence the realization of M. However, in many practical applications, H is 
nonergodic in one channel realization, and the AEM of HH* depend on the specific realization of M. 
This is the case to be investigated in this paper, where we assume that, although H is nonergodic in 
one channel realization, it is ergodic in the time domain. Specifically, we suppose that the realization 
of M is controlled by a set of random variables in 6; conditioned on a certain M, the AEM converge 
asymptotically to quantities that are functions of the realization of 0. We assume that H is ergodic in 
the time domain; thus, further averaging over the ensemble of 6 yields the unconditional AEM. 

The statistics of random matrix M are explicitly given below. We represent a random process by a 
function Z(u, t) of two variables, where u is a point that varies over the sample space U, and t is a point 
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over an index set T. In the current context, we suppose that a particular column k of M has a different 
statistical property from all of the others. This occurs when, for example, user k is the desired user in 
the processing of the linear model y = Hx + w. 

Definition 4: Consider a realization of random matrix M and let V(i,j) = |Mjj| 2 . For each given 
j G [1,K] \ {k}, {V(i,j) ■ 1 < i < N} is a sample function of the discrete-time stochastic process 
Zj(u,t) whose moments up to a certain order exist, and {Zj(u,t) : j G [1,K] \ {k}} is a collection of 
i.i.d. random processes. Furthermore, for the k-th column of M, {V(i, k) : 1 < i < N} is a realization 
of the random process Zk{u, t). 3 ■ 

Due to i.i.d. of Zj(u,ty$,, we will use statistics of Z(u,t) to represent those of Zj(u,t)\. By (1), for 
any y G {(j — 1)/K : 1 < j < K, j / k}, we can obtain a continuous-time random process, in which 
x < 1} is its realization. This continuous-time random process inherits properties from 
its discrete-time counterpart. That is, its moments up to a certain order exist, processes corresponding 
to distinct index y are i.i.d., and the statistics of these i.i.d. processes are characterized by the process 
z( N \u, t). When N — > oo, the limit of z^ N \u, t) is represented by z(u, t). A continuous-time counterpart 
z^ N \u,t) of Zk(u,t) can be found using (2), and the limit of z^\u,t) is symbolized by Zk(u,t) when 
N -> oo. 

We define three random sequences related to moments of the Wishart-type random matrix HH+: 



(JV) 
m 

AN) 



iV 1 tr{(HHt) m }, (5) 



€1 •■= m^n^, (6) 

C : = ht(H. fe Hl fe )-h fc , (7) 

where the superscript ^ denotes that the quantity is evaluated when H has N number of rows, is 
the k-th column of H, and is H with removed. Suppose that /jl^P, rj^l and 5^1 converge to 
l^m, Vm,k an d S mt k, respectively, when K, N — > oo and K/N — > (3. The importance of these asymptotic 
moment-related quantities are explained below. 

• The quantity fi m is the general definition for the m-th AEM of HH*. 



3 The relation between Zj(u,t), j / k, and Zk{u,t) may be stated below. Let the former and the latter be defined by a set 
of random variables and Ok, respectively. In many occasions, 6k C 0, and Zk(u, t) is obtained by fixing the random variables 
of 8 \ 0k by known deterministic values. The rationale of this lays on that the fc-th column corresponds to the desired user and 
thus more information is available. 
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• Given the model y = Hx + w described in Introduction, assume the k-th entry of x is of interest, 
and H is known to the receiver. Then r} m ^ is the quantity concerned with the detection of the desired 
symbol. 

• The quantity 5 m ^ is often used in the design and performance evaluation of a reduced-rank receiver 
with the received signal model y = Hx + w [10], [17]— [19], [33]. In [16], [34], i] m ^ is used for 
the same goal. 

In the following, some definitions regarding set partition theory are summarized. The properties that 
will be used in our later derivations are also provided. 

Definition 5 (Noncrossing Partition [35], [36]): Let S be a finite totally ordered set. 

1) We call w = {Bi, • • • , Bj} a partition of the set S if and only if Bi, 1 < i < j, are pairwise 
disjoint, non-empty subsets of S such that B\ U • • • U Bj = S. We call B\ , • • • ,Bj the blocks 
of w, and \w\ and \B{\ represent the number of blocks in w and the number of elements in Bi, 
respectively. The blocks B\ , ■ ■ • ,Bj are ordered according to the minimum element in each block. 
That is, the minimum element in B^ is smaller than that in B\ if k < I. For p G S, B m (p) denotes 
the index of the block that p belongs to under the partition w. For example, if p G Bi under w, 
then B m (p) = i. 

2) The collection of all partitions of S can be viewed as a partially ordered set (poset) in which 
the partitions are ordered by refinement: if w, a are two partitions of S, we have w < a if each 
block of w is contained in a block of a. For example, when S = {1,2,3,4,5,6,7}, we have 
{{1}, {2,5}, {3, 4}, {6}, {7}} < {{1,3, 4}, {2,5}, {6, 7}}. 

3) A partition w of the set S is called crossing if there exist pi < qi < P2 < 12 in S such that p\ 
and p2 belong to one block, and q\ and q2 belong to another. If w is not crossing, then it is called 
noncrossing. ■ 

The set of all noncrossing partitions of S is denoted by NC(S). In the special case of S = {1, ■ ■ ■ , m}, 
we denote this by NC(m). 

Definition 6 (Kreweras Complementation Map [35], [36]): Consider numbers 1, 2, • • • , m, and we form 
a totally ordered set {1, 1, 2, 2, •• • , m,m} by interlacing them with 1, 2, •• • , m. Let w G NC({1, ■ ■ ■ ,m}). 
Then its Kreweras complementation map K(zu) : NC({1, • • • , m}) — > NC({1, ■ ■ ■ ,fn}) is defined to be 
the largest element among those a G NC({1, • • • , m}) such that zuUa belongs to NC({1, 1, • • • , m, fn}), 
where "largest" is in the sense described in item 2) of Definition 5. ■ 

In this paper, noncrossing partition and Kreweras complementation map are employed to derive the 
limits of sequences (5)-(7). The same tool has been employed by [19], [37], [38] in obtaining AEM when 
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entries of H are i.i.d. A convenient representation of a noncrossing partition w as well as its Kreweras 
complementation map K(w) is the If -graph detailed in Appendix I. 



III. Asymptotic Moments Conditioned on a Certain Realization 
Consider an N x K random matrix H as described in Section II. We define 

,( ff )i„W/ T „n a.. , ■= U.W 



and 



£ m := ton WI/'^b)}, W== E{^>W(z,y)}, 

A,iV— >oo K,N—*oo 

K/N-+/3 



(8) 



In this section, we compute the above three quantities that are conditioned on a specific v^ N \x,y), and 
we show that fi m , f] m ^, and 5 m ^ are the limiting values of the sequences given in (5), (6) and (7), 
respectively, for each particular realization of v^ N \x,y). 

Theorem 1: Conditioned on v <yN \x,y), the random sequences Hm \ Vml an( * ^mfe converge a.s. to 
frm, Vm,k an d <5 m ,fc> respectively, as K, N — > oo and K/N — > /3. We have 



f M 

roeiVC(m) [ i=l ^efBsr (w) (t):t6B»} 

to={Bi,--- 



(9) 



f]m,k 



/3 |rohl E < 

weNC(m) 
ro={Bi,--- ,B|„|} 



H 



k Si6{Bjf (ro) (t):teBi} i=2s i 6{Bjf( ro) (t):t€B i } 



(10) 



and 



£ /?H E < 



roGA r C(m) 
ro ={B l5 -,B|„|} 



H 

II v(X Si ,Yi) 

«=1 s i 6{B K(w) (t):teB i } 



(11) 



where the expectations are with respect to continuous i.i.d. random variables X\, • • • , X m _^ + i, Yi, ■ ■ ■ ,Y\ Z 
uniformly distributed in [0, 1). 

Proof: The proof is provided in Appendix II, which relies on the representation of if-graph. ■ 
Note that, since the AEM of HH^ are independent of realizations of M if H is ergodic in one channel 
realization, the derived conditional AEM in Theorem 1 can also serve as the AEM formula for ergodic 
H when any realization of M is available. Issues regarding the complexity reduction in obtaining fi m , 
7/ m fc and 5 m ^k above are addressed before we move on to taking ensemble average of conditional AEM 
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according to the distribution of v^ N \x,y). To evaluate the expressions in (9)-(ll), it requires (m + 1)- 
dimensional integrals with the m-th Catalan number of times 4 . When the order m of the moment is high, 
the computation cost is huge. If the conditional variance profile v(x,y) possesses a special structure, the 
AEM expressions can be simplified through properties of noncrossing partitions. For instance, if v(x,y) 
is decomposable as g(x)h(y), where g(x) and h(y) are nonnegative valued functions with span [0, 1), 
then formulas given in Theorem 1 admit simpler forms. The complexity can also be reduced by resorting 
to recursions in conditional AEM expressions. We explore these two ways of complexity reduction in 
Theorems 2 and 3. 

The simpler form of fi m when v(x,y) = g(x)h(y) has been obtained in Theorem 3 of [10]. For 
completeness, formulas of fi m , fj m ^ and 5 m ^ given that v(x,y) = g(x)h(y) are all provided in the 
following theorem. 

Theorem 2: If the conditional variance profile is decomposable as v(x,y) = g{x)h{y) and Vk(x) = 
a k g(x), then 

1) fi m can be given by 




6i> - >6i>l ci>->c m _ !+ i>l 



i m-l+l 



x]jE{h(Xf} H E{g(Xy>}, (12) 




b^-^b^l Ci>->C m _ i + 1 >l 




(13) 



4 The number of noncrossing partitions of an m-element set is equal to the m-th Catalan number. 
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3) 5 mjk can be given by 

(m -/)!(/-!)! 



m, 

E E 



1=1 6l+ ... + 6 ; = m c 1 + ...+c m _ 1+1 =m^^'-- - ' 6 i)/(ci,"- ,C m _ m ) 
bi>->6i>l ci>->c m _ !+ i>l 
£ m— Z+l m— Z+l _j- ,„> c ^.j-, 

^nw) n e e/Sm • (14) 

i=l i=l n=l lyv ; 

Proof: See Appendix III. ■ 

In the following theorem, recursive formulas of fi m , fj mik and 5 mik are provided. They have simpler 
forms compared with those given in [10] for fi m and 8 m ^ k (Theorems 1 and 2 of [10] with notations 
fi m and <5 m> fc) by means of expanding the Stieltjes transform equations (3) and (4). For some particular 
variance profile, it is possible to obtain closed-form conditional AEM expressions by symbolic operations 
of scientific computation softwares such as MATLAB. 

Theorem 3: The a.s. limiting value of E{/j,^P\v^ N \x,y)} can be given as 

fi m = (3-E{il m (X,Y)}, (15) 

where the expectation is over i.i.d. random variables X and Y that are uniformly distributed in [0, 1), 
and pm{x,y) can be obtained from the recursion 

(m—l \ 
^ E E {Ai-i(^2/)} E {Am-i(a;,^)} + E{p m - 1 (X,y)}j , m>l (16) 

with jj,o(x,y) = 1. Also, we have 

f] nhk = E{fj mjk (X)} and 5 m , k = /? • E{v k (X)jl m (X, Y)}, (17) 

where X and F are i.i.d. uniform random variables in [0, 1), and rj m ,k{x) can be obtained recursively 
from 

/ m-1 \ 

Vm,k( x ) = v k (x) (3 ^2 Vi-i,kE{iji m -i(x, Y)} + J7 m _i jfc , m > 1 (18) 



i=i 



with 77o,fc(^) = 1- 

Proof: See Appendix IV. 



IV. Unconditional AEM 



Unconditional AEM can be obtained as the ensemble averages of conditional AEM fi m , fj mjk and 
<5 m fc according to the distribution of v(x,y). It is readily seen that performing the above operation to 
recursion based AEM expressions is inconvenient. Thus, it is suggested that results in Theorems 1 and 2 
are employed. We denote the ensemble averages of fi m , fj mjk and 5 m ^ k by \i m , rj m:k , and S m:k , respectively. 
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Recall, in Definition 4, we describe that, for any j £ \ {k}, {V(i,j) : 1 < i < N} is a 

realization of random process Zj(u,t), and Zj(u, t)'s are i.i.d. for distinct j. We use the statistics of 
Z(u,t) to denote those of Zj(u,ty&. If Xo is a set of deterministic integers, we define Momz(Xo) as 
the |-X*o|-th moment of Z(u,t) at Xq, i.e. 



Mom z (X ) =e| 11 Z(u,x) I. 



(19) 



Moreover, if X = {Xl,-- - , X n } is a set consisting of random variables and X C X, we define 
Mom^|x (Xq) as the moment of Z(u,t) conditioned on X, i.e. 



X\ , • • • , X n 



(20) 



Mom Z | X (X ) = H Z(u,x) 

The definitions of moments in (19) and (20) hold also for random processes Z^{u, t), (u, t), (u, t), 
z(u,t) and Zk(u,t) depicted in Section II. 

Theorem 4: The unconditional AEM Hm \ r^\, and converge a.s. to 



Mr, 



f]m,k 



w€NC(m) 
H7={S 1 ,---,_B|„|} 



M 
i=l 



roGA r C(m) 



H 



(21) 



j=2 



M 



i=i 



roGAfC(m) 

respectively, where X = {X±, ■ ■ ■ , X m _| ro | +1 } is a set of i.i.d. random variables uniformly distributed 
in [0, 1), and = {X l : I G {£^ M (t) : t G B t }}. 

Proof: See Appendix V. ■ 
Two examples of nonergodic random channel matrices are considered in the following subsections. 



A. Variance Profile with Each Column a Switching Process 

Let H be a random matrix as described in Section II. Suppose that elements at the first column of 
H have the same variance (nonzero in general); for any other column j / 1, {V(i,j) : 1 < i < N} 
is a switching function from to a positive value or vice versa with a random switching time, where 
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the positive value is random. This random matrix H corresponds to the channel matrix of a one-shot 
asynchronous CDMA to be detailed in Section V-B. 
We describe V(i, j) as follows. For j = 1, 

V(i,l)=m(l), l<i<N. (22) 

For j 6 [2, if], {V(i,j) : 1 < i < N} is a switching function from to m{j) or vice versa, given by 

\rc -\ \ m (3>{i-r{j)), if w(j) = 0, 

= < 1 < i < AT, (23) 

[ m(j)[l - u(i - r(j))], if w(j) = 1, 
where m(j) is positive and governs the magnitude, u(i) is the unit step function, r(j) € [0,N — 1] 
specifying the switching time, and w{j) G {0, 1} controls the states before and after switching. Suppose 
that {m(j) : 1 < j < K}, {r(j) : 1 < j < K} and {w (j) : 1 < j < K} are realizations of independent 
random variables M, T and W, respectively. In particular, T is uniform, and W is equal to or 1 
equally probably. It is clear to see the asymptotic empirical distribution of V(l,j), V(2,j),--- , V(N,j) 
depends on the realizations m(j), r(j) and w (j); thus, the random matrix H is nonergodic in one channel 
realization. 

Theorem 5: Given the random channel matrix described in (22) and (23), the asymptotic m-th moment 
fi m is expressed as 

zueNC(m) i=l 
ro={Bi,-,B M } 

the asymptotic m-th moment i] m ^ is given as 



M 

Q (1 - maxX^-j + min X m -,j) 
j'=i 



(24) 



M (IH 

w= e (/3/2) H - i n E { M ' Bi 7 E i n^-^^+^^i) f ' (25) 

weNC(m) i=l \j=2 ' 



and 



where is as defined in Theorem 4. 
Proof: See Appendix VI. 



<5 m ,i = E{M}^ m , (26) 



6. Variance Profile with Each Column a Bernoulli Process 

Consider a random matrix H as described in Section II. Let the elements at the first column of H have 
the same variance (nonzero in general); for any other column j / 1, {V(i,j)} is a Bernoulli process 
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having states m(j) > and with probabilities p(j) and 1 —p(j), respectively. In specific, when j = 1, 
V(i, 1) = m(l); for each j e [2, K\, V(i,j) for 1 < % < N are i.i.d. random variables having the density 

fj(v) = {l-p(j))5(v)+p(j)S{v-m(j)), 0<p(j)<l, m(j)>0, 

where {m(j) : 1 < j < K} and {p(j) : 2 < j < K} are realizations of independent random variables 
M and P, respectively. 

Theorem 6: For the variance profile composed of Bernoulli processes as defined above, the asymptotic 
m-fh moment fi m is equal to 

<- = E^ E ^£^^fl^' h M^}. 07) 

1=1 b 1 +b 2 +--+b l =m ' * 7 i=l 

&i>6 2 >->&i>1 

?7 m) i is equal to 

Z=l 6l+ ... +6i=m ' '/ i=l n=l 1 J 

6!>--->6 ; >l 

and 

<5 m ,i = E{M}^ m . 

Proof: See Appendix VI. ■ 
In the following, the relation of AEM's yielded by variance profiles of switching processes and a 
Bernoulli process is built by some approximations. We use formulas of \i m in (24) and (27) as an 
example. When the expectation 

E {n ( l- max X^,j + min X w .j) J (29) 

in (24) is approximated by 

M M / | B | 1 n M 2 

[]E{1 - max^j + min^ ro; ,} = ]J I 1 - -J-^ + r^TT-j J = II j^ppy, (30) 

j=l j=l ^ / j=i 

we obtain (24) approximately as 

*. * £>' e m( "'m; ) fe ;. ( " 1 M +2> ri E{Mtil( ''- +irl - (31) 

1=1 b 1 +b 2 +-+b l =m A°l>°2> ,Oj; i=1 

where the first equality of (30) holds because Xj's are i.i.d. uniform random variables. Note that, when 
all but one of BiS are singletons, (29) and (30) are equal. 
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Consider the special case that, in the variance profile of Bernoulli processes, the random variable P 
is uniformly distributed in (0, 1). Then, E{P n } = (n + l) -1 , and the asymptotic moment /x m in (27) 
becomes the same as (31). The coincidence of the two expressions can be interpreted as follows. For 
the variance profile as a switching process, each column of H has a continuous segment of zeros with 
its length uniformly distributed in [0, N). If these zeros are spread out randomly to the whole interval 
[0, N), the variance profile of switching processes becomes a variance profile of Bernoulli processes. 



In this section, two applications of the AEM formulas derived above are provided. In the first subsection, 
the spectral efficiency of an MIMO channel given a conditional variance profile is computed. Our intention 
is to exemplify the use of the moment based method, instead of the widespread Stieltjes transform based 
approach, in determining the spectral efficiency of a communication system. In the second subsection, 
we consider a one-shot symbol-asynchronous but chip-synchronous CDMA system, where the random 
variance profile depicted in Section IV-A (as well as IV-B, approximately) is the corresponding variance 
profile. 

A. Correlated Fading in MIMO 

In some applications of probability, it is frequent that the (infinite) moment sequence of an unknown 
distribution F is available, and these moments determine a unique distribution. Suppose that the final aim 
is to calculate the expected value of function g(X) of the random variable X whose distribution F is 
unknown. One of the most widely used techniques for evaluating E{g(X)} is the Gauss quadrature method 
[15], where 2Q + 1 moments {m„}^L of distribution F are used to determine a Q-point quadrature rule 
{w q ,x q }® =1 such that 



and the approximation error becomes negligible when Q is large. In this subsection, the Gauss quadrature 
method is used to compute the mutual information of a spatially correlated MIMO system. 

Consider re T transmit and ur receive antennas with the corresponding multiantenna channel denoted 
by H, whose (i,j)-th entry is the fading coefficient between the j-th transmit antenna and i-th receive 
antenna. Let the correlation between the (i, j)-th and (Y, j')-th entries of H be represented by 



V. Applications 




(32) 



Rn(i,j;i',j 



■i 



) - E{HijH*, tj ,}. 
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It is shown in [3] that any MIMO channel with correlation i?H that falls within the unitary-independent- 
unitary (UIU) 5 framework and has bounded eigenvalues is unitarily equivalent to an i.n.d. channel H with 
the mean and variance of the (i,j)-th entry equal to and the (i,j)-th eigenvalue of i?H, respectively. 
That is, HH f and HH* have the same AED. 

We consider Example 2 of [3], where n T = 3, n R = 2, and an i.n.d. channel H whose component is 
zero-mean and has variance equal to the entry of G at the same location with 

0.4 3.6 0.5 
0.3 1 0.2 

each (i, j)-th entry of H has independent real and imaginary parts uniformly distributed in the interval 
[—y/1.5Gi t j, y/1.5G{j]. It is known that the normalized input-output mutual information is 



G = 



(33) 



— ^log 2 (l + 7 A, (H*Ht))^> J log 2 (l+ 7 X)dF(X), 

where 7 is the input SNR, Aj(-) is the i-th eigenvalue of the indicated matrix, X is the limiting random 
variable governing the eigenvalue distribution of H3>ff with <1> the diagonal input covariance matrix, 
and F(X) is the cumulative distribution function of X. Thus, the liming mutual information can be 
obtained by letting g(x) = log 2 (l + jx) in (32), and the AEM of H«&Ht can be acquired from those of 
HHt and * using the result of [39]. 

Fig. 1 shows both the simulated and analytical results of the capacity and the spectral efficiency with 
uniform power allocation versus input SNR, where the capacity-achieving power allocation $ is found 
using the algorithm in [40], and two curves showing analytical results (solid and dashed lines) are obtained 
using a 5-point Gauss quadrature rule. The analytical results in Fig. 1 using the moment-based method 
have a slight performance gain compared with the transform-based method presented in [3]. The gain in 
terms of error percentage 6 becomes more and more obvious when the input SNR gets larger. When the 
input SNR is equal to 20dB, capacities based on the Monte-Carlo simulation, the methods in [3] and this 
paper are 11.49, 11.09, and 11.41 (bits/s/Hz), respectively; spectral efficiency with uniform power of the 
three methods in the same order are 11.37, 10.99, and 11.17 (bits/s/Hz). 



5 For the definition of UIU, see [3] 

6 The error percentage is defined as (simulation result-analytical result)/simulation result x 100%. 
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Fig. 1. The capacity and the spectral efficiency with uniform power allocation of a spatially correlated MIMO with variance 
profile given in (33). The two analytical curves (solid and dashed lines) are obtained by 5-point Gauss quadrature method. 



B. One-Shot Asynchronous DS-CDMA 

Consider a DS-CDMA system with asynchronous transmission where each user's spreading sequence 
is chosen randomly and independently. It is assumed that the system is chip-synchronous, i.e., the relative 
delay of each user is an integer multiple of the chip duration. 

Suppose that b k £ {+1, —1} is the transmitted symbol of user k, s fc g {+1/^,-1/^}^ is the 
signature sequence of the same user, and n is the AWGN with variance a 2 . We adopt the signal model 
analogous to that presented in [8, eqn. (2)]; the sampled discrete-time model for the received signal is 
given as 

K K 

r = aiMi + ^2 akh^k + ^ a k d k N k + n, (34) 

k=2 k=2 

where a k is the fading coefficient of the fc-th user, b k ,d k G R are two consecutive symbols of the k-ih 
user which overlap with user 1 in the observation window. Symbols b k and d k have effective signature 
sequences u^. and v&, respectively. If t k £ Z + denotes the relative delay of user k to user 1 in terms of 
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(sk) N -t k +v 1 <i< t k 

and (vfc)i = < 

0, t k <i<N 



0, 1 < i < t fc 
(sfc)i_ tfc , t k <i<N 



the number of chips, then 

The signal r in (34) can be expressed in a more compact form as 

r = SAb + n, (35) 

where the notations are defined as follows: 

S = [ui, u 2 , • • • , u Kl v 2l ■ ■ ■ , v K ] € R^ 2 *" 1 ), 

A = diag(a 1 ,« 2 ,-- - ,a K ,a 2 ,--- .a^eC^W 2 ^ 1 ', 

b=[6i,62,--- ,b K ,d 2 ,--- ,4] T Gl (2if - 1)xl . 

Define H = SA and R = H^H. Following the methodology of [16], [34], we obtain a D-dimensional 
reduced-rank MMSE receiver for user 1, i.e. the estimate of b\ as &i(r) = w T r, where 

w = 3>~ V, (36) 

with 

<p = [(R)n, (R 2 ) n , • • • , (R D )n] T and = (R+ J )n + er 2 (R i+J ' -1 ) 11 , l<i,j<D. (37) 

It is seen (R m )n is equal to r^-j* defined in (6). 

Suppose that the relative delays t 2 , • • • ,tx are realizations of a discrete random variable uniformly 
distributed in [0, iV). The random matrix H = SA has a conditional variance profile identical to that 
presented in Section IV- A, where |«i| 2 ,--- ,|ok-| 2 are realizations of random variable M. The AEM 
formula rj mj i given in (25) can be used to obtain the reduced-rank MMSE receiver w = Q l tp in 
the asymptotic regime. However, the computation of r\ m ^\ is involved with several multi-dimensional 
integrals, and the complexity is quite high. We may resort to its counterpart in (28) whose variance 
profile is composed of Bernoulli processes, and the random variable P therein is set as uniform in (0, 1). 
As demonstrated in the discussion following Theorem 6, the use of (28) as a replacement of (25) is 
an approximation by spreading the zero elements in each column of channel matrix H randomly and 
uniformly. 

Simulations are run to compare the AEM's by (25) and (28). The AEM of the former are obtained 
through 10 5 -point Monte-Carlo integration, and we evaluate the first eight moments for various values 
of (5. Since the distinction of the two sets of AEM is on the distribution of zero entries in columns of 
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2 4 6 8 10 12 14 16 18 20 

Input SNR (dB) 



Fig. 2. The output SINR of three receivers in a symbol-asynchronous but chip-synchronous CDMA system with an unfaded 
channel: K = 16, N = 64, and D = 4. 

H and irrelevant to the fading process M, we let M as a constant equal to 1. The random variable P 
in (28) is set uniform in (0, 1). The results show the error percentage 7 increases with the order m of 
moments, equal to 0.02%, 0.22%, 0.90% and 1.78% for m = 2, 4, 6 and 8, respectively. It is seen that 
AEM obtained from (25) have larger values than those from (28). 

In the following, we compare the output SINR of three receivers: i) the full-rank MMSE receiver, 
ii) the reduced-rank MMSE receiver that uses the actual weights, and iii) the asymptotic reduced-rank 
MMSE receiver that uses the asymptotic weights. For the second receiver, the filter in (36) is employed; 
while for the third receiver, entries of $ and (p in (37) are replaced with their limiting values in the 
asymptotic regime. The number of users is K = 16, and the spreading gain is N = 64. The two receivers 
have the rank D = 4. For the asymptotic weights of receiver iii), we use (28) to compute r/ m i because 
it yields a very close result to (25) and it demands lower complexity. The output SINR is the ratio of 
total signal power and total interference -plus-noise power of a large number of independent simulations. 



7 The error percentage is defined as — — — — x 100%. 
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2 4 6 8 10 12 14 16 18 20 

Input SNR (dB) 



Fig. 3. The output SINR of three receivers in a symbol-asynchronous but chip-synchronous CDMA system with a Rayleigh 
flat fading channel: K = 16, N = 64, and D — 4. 

Fig. 2 shows the output SINR for the three receivers in an unfaded channel. It is seen that the output 
SINR of the asymptotic reduced-rank receiver is close to that of the receiver using actual weights for 
the input SNR considered, although the asymptotic weights of the third receiver are obtained using the 
ensemble statistics of relative delays and spreading codes. 

Fig. 3 plots the performance of the three receivers in a Rayleigh flat-fading channel. In formula (28), 
the moments of M are obtained as those of a complex Gaussian with independent real and imaginary 
parts. As the empirical distribution of the fading coefficients cci,- ■ ■ ,<xk can hardly approach Gaussian 
in one channel realization, the asymptotic weights can be viewed as calculated using a fading statistic of a 
very long observation time. It is seen from the figure that the output SINR of the asymptotic reduced-rank 
receiver is larger than the other two receivers until the input SNR exceeds 16dB. This is because, for the 
asymptotic receiver, fading coefficients of extremely large magnitudes that occur occasionally result in 
profound effects on the output signal power. For the other two receivers, the effect of fading coefficients 
with large magnitudes can be offset by the MMSE filtering that is computed from the known fading 
coefficients, while the asymptotic receiver fails to do so because its weights are computed based on the 
ensemble fading process over a long time. In Fig. 4, we show the bit error rates (BER) versus input SNR 
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Fig. 4. The BER of three receivers in a symbol-asynchronous but chip- synchronous CDMA system with a Rayleigh flat fading 
channel: K = 16, N = 64, and D = 4. 



of the three receivers in the Rayleigh fading environments. Although the asymptotic reduced-rank receiver 
has the largest output SINR, its BER performance is worse than the other two receivers as expected. To 
boost the performance of the asymptotic receiver, we can employ the implementation in [10] that the 
moments of the fading process are found by the empirical distribution of finite number of realizations. 

VI. Conclusions 

In this paper, we have used noncrossing partition to derive expressions of several AEM related quantities 
of a large Wishart-type random matrix HH* when the variance profile of H is random and H is 
nonergodic in one channel realization. These quantities are useful in the design and analysis of a number 
of communication systems. It was assumed that H = M S, where entries of S are zero-mean i.i.d. 
random variables, and columns {|Mjj| 2 : 1 < % < N} are realizations of i.i.d random processes. The 
derivation started from obtaining the conditional AEM given a certain realization of M, and we then 
derived unconditional AEM by computing ensemble average of the conditional quantity according to the 
distribution of M. We have obtained conditional AEM expressions in both direct and recursive forms. 
When the conditional variance profile is decomposable as the product of two one-dimensional functions 
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with deterministic nonnegative values, the conditional AEM expressions can be simplified. Two statistical 
models of M were given. One is that each column {|Mj j\ 2 : 1 < i < N} is a switching function with 
a random switching time; the other is that {\M{ j\ 2 : 1 < i < N} is a Bernoulli process with a random 
success probability. 

Two application cases have been provided. One is the use of Gauss quadrature method to computing 
the ergodic capacity of a spatially correlated MIMO system by means of AEM; a better accuracy than 
the Stieltjes transform based method of [3] has been achieved. The other is the design of a reduced-rank 
MMSE receiver in a one-shot asynchronous CDMA system, whose corresponding variance profile has 
each column as a switching function with random switching time; the AEM calculation is involved with 
multidimensional integrals. Due to the difficulty in evaluating the integrals, a variance profile consisting 
of Bernoulli processes was used instead. Numerical results showed the approximation error is negligible. 

Appendix I 

if-GRAPH 

Definition 7 (K -graph [41]): The .KT-graph of a noncrossing partition w = {B\, ■ ■ ■ , Bj} G NC(m) 
is denoted by G = (V,£). The vertex set is V = {1,2, • • • , j}, and the edge set is £ = {1,2, ■ ■ ■ ,fn}. 
Edge f connects vertices s and t if £> ro (r) = s and £> ro (r + 1) = t (with m + 1 := 1). ■ 
The if-graph ofro = {{l,5,7},{2,3,4},{6}} € NC(7) is shown in Fig. 5. The properties of a K-graph 
are summarized below. 

1) A noncrossing partition can be recovered from its K-graph. Blocks of w can be identified by 
starting at vertex 1 and traversing edges 1, 2, • • • ,m; if edge r starts at vertex i, then r G B; L . For 
example, in Fig. 5, edges (1, 2, • • • , 7) start at vertices (1, 2, 2, 2, 1, 3, 1), respectively. 

2) Due to the noncrossing property, a i^-graph of w G NC(m) and \m\ = j is a concatenation of 
m—j + 1 cycles with any pair of cycles connecting by at most one vertex. The K-graph in Fig. 5 
is a concatenation of 7 — 3 + 1 = 5 cycles, denoted by C±, ■ ■ ■ , C5. The cycles are ordered in the 
ascending order of the minimum edge element in them. 

3) The z-th block of K(w) is recognized by the edges in the i-th cycle of the i"C-graph. In Fig. 5, 
since edges {1,4}, {2}, {3}, {5,6}, and {7} constitute cycles C\ , C2 , C3 , C4 and C5, respectively, 
we have K{w) = {{T, 4}, {2}, {3}, {5, 6}, {7}}. 

4) There is a bijective correspondence between blocks of a noncrossing partition w and vertices of 
its X-graph. There is a bijective correspondence between blocks of the Kreweras complementation 
map K(w) and cycles of the K-graph of w. 
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Fig. 5. The A"-graph of noncrossing partition vj = {{1, 5, 7}, {2, 3, 4}, {6}}, which is composed of five cycles Ci, • • • , C5. 



Remark: A if-graph can be interpreted in a more visually convenient way as follows. Let us arrange 
vertices v\ , ■ ■ ■ ,v m orderly (either clockwise or counter-clockwise) in an m-vertex cycle, and let edge r 
connect vertices v r and tv+i- The K-graph of w € NC(m) can be obtained by, for 1 < i < \w\, merging 
vertices {vj : j E Bi} into a vertex i. When v/s are merged into a new one, edges originally incident on 
vfs become incident on the new vertex. Mergence of two adjacent vertices results in a self-loop cycle. 
■ 

Take Fig. 5 as an example. Vertex 1 is formed by the mergence of vertices v\, v§ and vj of a 7- vertex 
cycle. Edges 1 and 7, 4 and 5, and 6 and 7 originally incident on vertices v%, v§ and v-j, respectively, 
become incident on the new vertex 1 after mergence. Moreover, since vertices v\ and vj are adjacent in 
a 7-vertex cycle, their mergence results in a self-loop cycle C5. 

Appendix II 
Proof of Theorem 1 

A. Proof for fi m 

For the ease of extending the derivation of fi m to f) m ^ and b m ^, we rewrite the right-hand-side of (5) 

as 

f3 ■ K~ x tr{(H^"H) m }, (38) 
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(a) (b) 



Fig. 6. The /("-graphs of (a) m-block noncrossing partition {{1}, {2}, • ■ ■ , {m}}, and (b) (m — l)-block noncrossing partition 
{{1}, {2}, ■■• ,{s, t}, {s + 1}, ■ • • {m}} with s < t. 

and the operand of the limit in the expression of fi m is given as _1 E{tr{(H"'"H) m }|?/ iV )(x, y)}. For 
simplicity, the conditional notation {-\v ( N ) (x, y)} is omitted below. We have 

K _1 E{tr{(H^H) m }} 

= K- 1 E{(HtH) fclfc2 (HtH) fc2fc3 ---(HtH) fcmfcl } 

{fei,-,fc m }e[i,A'] m 

= K 1 ^2 E {( H n m M Hr >-uki) { H n u k 2 H nv,k 2 ) " " " [ H n m - 1 ,k m H n m ,k m ^ ■ (39) 

{k x ,- ,k m }e[l,K}™ 
{m,- ,n m }e[l,N]"> 

To evaluate the expectation in (39), it is required to consider the equivalence relation of variables 
k\, ■ ■ ■ , k m as well as that of n\, ■ ■ ■ , n m . As the equivalence relation is equivalent to set partitioning by 
that two variables take the same value ( in the range of [1, K] or [1, N] depending on whether ki 's or rii 's 
are considered) if and only if they are partitioned in the same block, set partition theory can be employed 
as the tool for expectation evaluation of (39). Thus, we investigate partitions of two totally ordered set 
{fci, • • • , k m } and {ni, • • • , n m }. It is known that, when evaluating (39) asymptotically, it is sufficient 
to consider the situations that k\ , • • • , k m is partitioned noncrossingly. We will show that, to yield non- 
vanishing expectation of (39), the partition of {m, • • • ,n m } should be the Kreweras complementation 
map of the partition of {&!,••• ,k rn }. In the following, for convenience, we denote the partitions of 
{ki, k2, ■ ■ ■ , k m } and {m, 712, • • • , n m } by those of {1, 2, • • • , m} and {1, 2, • • • , m}, respectively. 

We consider w = {B\, • • • , Bj} G NC({1, • • • , m}) for various values of j. First, we consider j = m, 
resulting in w = {{1},{2}, • • • ,{m}}. By the connection of equivalence relation and set partitioning, 
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this corresponds to k\ ^ ■ ■ ■ ^ k m in (39). In this case, since entries of H are independent, (39) becomes 

K- 1 E {K m M H niMMK lM H n 2 ,k 2 } ■ ■■m* nm . 1 ,k m Hn m ,k m }. (40) 

all ki's distinct 

Til )*** 1TI771 

To yield a nonzero summand in (40), it is required that n\ = ■ ■ ■ = n m . Interpreting this equivalence 
relation by partition of {m, • • • ,n m }, we can see it corresponds to the Kreweras complementation map 
of {{1}, • • • , {w,}}, i.e. if ({{1}, • • • , {m}}) = {{1, • • • ,m}}. The sum in (40) thus becomes 

E E fl^ («) 

{b u - ,b m }€{l,K}™ Cl e{l,N] i=l 

The If -graph of {{1}, • • • , {m}} is shown in Fig. 6(a), where there are m vertices 1, • • • , m, and the 
if -graph is itself a cycle labelled as C±. The variable bi in (41) is interpreted as the integer in [l,if] 
chosen by the z-th block of the partition of {{1}, • • • , {m}} (or vertex i of the if -graph), and c\ in (41) 
is the integer chosen by the block of {{1, • • • ,m}} (or the cycle of the if -graph). In the sequel, bi (or 
Cj) is used to represent the integer in [l,if] (or [l,iV]) selected by vertex i (or cycle Cj). 

Next, consider w = {B\, • • • , Bj} G NC({1, ■ ■ ■ , m}) with j = m — 1, where one of the blocks 
contains two elements s and t (s < t), and all other blocks are singletons. In this case, (39) becomes 

K ^ ^2 E { i ^n m ,fe 1 - ff n 1 ,fe 1 }E{ii'* life2 ii'n 2 ,fe 2 }- ■ ■ 

all ki's distinct except for k s = k t 

,n m 

E { H n s - U k s H n s ,k s Hn t _^k t H n t ,k t } " " " E { H n m - 1 ,k m H n m ,k m } ■ (42) 

To yield a nonzero summand of (42), we need that the variables n\ , ■ ■ ■ , n m in each expectation of (42) 
are paired. That is, n m = n.\ in the first expectation, ri\ = ri2 in the second, and so on. In the expectation 
containing n s _i, n s , nt-i and n t , there are three possible cases for pairing, i.e., (i) n s _i = n s and nt-i = 
n t , (ii) n s = n t and n s _i = nt-i, and (hi) n s _i = n t and n t -\ = n s . Among these, case (iii) (along 
with n m = n\, n\ = ri2, and so on) divides m, ■ ■ ■ , n m into two groups (m, • • • , n s _i, n<, • • • , n m ) and 
(n s , ■ ■ ■ , nt-i) with n/s in the same group take the same value. On the contrary, cases (i) and (ii) result 
in only one group. Thus, when computing (42) in the limit N — ► 00, it is sufficient to consider only case 
(iii) since it yields the highest dimension of N. Note that the grouping of rij's resulting from case (iii) can 



be obtained from if ({{1}, {2}, • • • , {s, t}, ■ ■ ■ {m}}) = {{!,■■■ , s - l,t, ■ ■ ■ ,m}, {a, ■ ■ ■ ,t - 1}}. The 
if -graph corresponding to the current case is shown in Fig. 6(b). We notice that case (iii) is equivalent 
to letting edges in the same cycle be in the same group. By assigning bi and Cj to the i-th vertex and 
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the j'-th cycle, respectively, (42) is given by 

K -i sr sr v< < c MV{ciM) v{c u b s )v{c 2 ,b s ) v{ Cl ,b m -i) 

2^ 2^ N at'" at2 "•• N ' ^> 

bi,---,6 m _i Ci,c 2 

where V(ci, b\)/N, V(ci,b2)/N and V(ci,b m -i)/N are equal to the first, second and fourth expectations 
of (42), respectively, while V(c\, b s )V(c2, b s )/N 2 comes from the third expectation. The expression of 
(43) can be recovered directly from Fig. 6(b). In specific, each vertex % is associated with 

n V(cj,bj) 
N ' 

ieT t 

where j E Ti if vertex i resides in cycle Cj, and (43) is obtained by multiplying the terms corresponding 
to each vertex altogether. For example, in Fig. 6(b), vertex s resides in cycles C\ and Ci, while vertex 
1 resides only in C\. Thus, T\ = {1} and T s = {1, 2}. We can obtain (41) from Fig. 6(a) in the same 
manner. 

In the above, we consider noncrossing partitions with m and m — 1 blocks. Here we summarize the 
findings in the two cases and extend to noncrossing partitions with arbitrary number of blocks. According 
to the remark following Definition 7, the if-graph of any noncrossing partition w of {1,- • • , m} can 
be obtained by successively merging two vertices of a same cycle. 8 This vertex mergence starts from an 
m-vertex cycle (Fig. 6(a)) where all edges are assigned with an identical integer in [l,iV]. At the first 
iteration, two vertices of the m-vertex cycle are merged into one, yielding a concatenation of two cycles 
(Fig. 6(b)). Edges in this two-cycle graph take the same integer if and only if they are in the same cycle. 
At the r-th iteration, 1 < r < m — \w\, two vertices s and t in any one of the r cycles are merged to 
yield r + 1 cycles in total. The cycle in which vertices s and t originally locate is then torn into two, 
and the edges in these two cycles are assigned with the same integer if and only if they are in the same 
cycle. 9 This assignment of rij's leads to that partitions of {k±, • • • , k m } and {ni, • • • ,n m } in (39) are 
Kreweras complementation map of each other. The contribution of w to (39) can be obtained from its 
K-graph as well as the integer assignments of fcj's and n^'s, given by 

K -i N -m £ n v ^m)--- n v( Cs ^,b lM ), (44) 

6i,---,6|w| ci,— ,c m _|„|+i si£Ti S| w |6T| ro 



If two vertices of distinct cycles are merged, it results in a crossing partition. 

9 The reason is stated right after (42), i.e., it has the highest dimension of N in the combinatorics of m's. 
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where Tj, 1 < i < is the set composed of the indices of cycles that vertex i resides. For example, for 
the noncrossing partition associated with Fig. 5, we have T\ = {1,4, 5}, T2 = {1, 2, 3} and T3 = {4}. 

We have the following key observation. Suppose that the i-th block of w contains k elements and 
is given by Bi = {Bi(l), • • • , Bi(k)}. Then, all the edges • • • , £>«(£;)} in the if-graph must be 

located in distinct cycles, and these k cycles are exactly the cycles that vertex i resides. For instance, 
the noncrossing partition in Fig. 5 has B\ = {1,5,7}. It is seen that edges 1,5 and 7 are in distinct 
cycles Ci,C4 and C5, respectively, and these three cycles are the ones that vertex B\ is located. Since 
the edges partitioning of {1, • • • ,fn} is the Kreweras complementation map of w for {1, • • • ,m}, the 
index of the cycle that edge r is located can be identified by Bk(zv) ( r ) (defined in Definition 5). Thus, 
given Bi = {Bi(l), ■ ■ ■ , Bi(\Bi\)}, the set T« in (44) can be expressed as 

T t = {B KM {Bi(l)),.-- ,^ M (^(|A|)} 

= {B K {vo){t) ■ t G Bi}, l<i<|o7|. 

It follows that (44) can be rewritten as (45) shown in the following lemma. 

Lemma 1: Suppose that w = {Bi,-- - ,-B| ro |} is a noncrossing partition of the totally ordered set 
{k\, • • • , k m }. Let fcj's in (39) take the same integer in [1, K] if and only if they are partitioned in the 
same block of w. Then, the sum of (39) with the highest order of N is given by 

M 

K -i N -m £ e n n (45) 

61,- ,b\„\e[l,K] ci,- ,c m _|„|+ie[l,JV] i=l s i e{Bjf(w ) (t):t6S i } 

bl^ — ^6|w| Cl^ — ^Cm-|w| + l 

Proof: The proof is straightforward from the preceding discussion. ■ 
We evaluate (45) asymptotically. Rewrite (45) as 

ci,-,c m _|„|+i i=l 6, s»6{Sjf( ro )(t):teBi} 

where the term inside the bracket {•} is identified as the Riemann sum of 

n « m 

Sie{B K ^)(t):teBi} 

in yi e [0, 1] with the partition and tag both being {i/K : < i < K — 1}. Similarly, the second line of 
(46) is the multi-dimensional Riemann sum of 

n/ n v {N \x Si ,y r )d yi 

i=l J ° s 4 €{S jr( „ ) (t):t6B i } 
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on {xi, ■ ■ ■ ,x m _| ro | +1 } G [0, l] m M+ 1 with the partition and tag on each dimension {i/N : < i < 
N — 1}. Evaluating (46) asymptotically, i.e., K, N — > oo and If/iV — ► /3, we obtain 

0M-1 /" ■■■[ <H[ H v(x ai ,y i )dy i \dx 1 dx 2 ---dx rn _ H+1 . (47) 

To get /i m , we sum up (47) for all G NC{m) and multiply the sum by (3 (due to (38)). By Theorem 1.1 
of [42] (or Lemma 2.60 of [12] in brief), we know [iffl converges a.s. to jj, m . 

B. Proofs for f)„ hk and <5 m , fc 

To show that fj mjk is equal to (10), we set variable b\ in (45) equal to k and multiply (45) by K. 
By writing the resultant equation in the form of Riemann sum and taking the limit K, N — ► oo and 
K/N -> /3, we can obtain (10). 

Next, we show that 5 m ^ k is equal to (11). The operand of the limit in the expression of 6 m:k is given 

as 

E { H n 1 ,k H n u k 2 H* 2! k 2 H n2jk3 - ■ ■ H* m+ukm+i H nm+uk }. (48) 

k 2 ,~ ,k m+1 e{l,K]\{k} 
tii,- ,n m+ i6[l,JV] 

Since &2, ■ ■ ■ , k m+ \ are unequal to k, we should let n\ = n m+ \, otherwise, the summand is zero. Thus, 
(48) can be further written as 



N- 1 



£ V(n lt k)E { (H nuk2 H* n2M ) ■ ■ ■ (H nm , km+1 H* nukm+i ) } . (49) 

k 2 ,- ,k m+1 e[i,K]\{k} 

Tii,"- ,n m 

Let us change the variable k m+ \ as fci. Then (49) becomes 

N- 1 £ nni,k)E{{H nmM H* niM ) {H niM H: 2M )...{H nm _ ukm H* nmik J}, (50) 

fci,-,fc m e[i,jq\W 

ni,— ,n m 

where the term inside expectation is the same as that of (39) except for a conjugate operation. Using a 

similar reasoning as we have proceeded to obtain (45), we can write (50) as 

H 

N -(m + l) £ V( Cl ,k)H H V^A), (51) 

bi,-,b\„\ ci,-,c m _| w |+i »=ls 4 €{Bjr(w)(t):*6Bi} 

&l^-^&|w| Ci^-^C m _| ro | + i 

with 6i, ■ • • , unequal to A;. Writing (51) as a Riemann sum, letting N — > oo, and summing up for 
all cc7 G NC(m), we can obtain (11). 

The proof of the a.s. convergence of random sequences rj^l and 5^1 follows from Theorem 1.1 of 
[42]. 
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Appendix III 
Proof of Theorem 2 

Part 1) of this theorem has been proved in [10]. Since the proofs of Parts 2) and 3) are very similar, 
here we only prove Part 2). The proof of Part 3) simply follows the same line. 

Consider w G NC(m) such that w = {Bi, • • • , B{\ and K{xu) = {C\, • • • , C m _/ + i}. When v(x, y) = 
g{x)h{y) and vu{x) = dkg{x), due to the property of noncrossing partition, the contribution of w to 
(10) can be written as 

l m-l+l 

P^a^HElHY^} [J H9(Xj) lCA }- (52) 

i=2 j=l 

Define a permutation operator V that yields a nonascending order sequence. For all w = {B±, ■ • • , B{\ G 
NC(m) such that V(\B\\,--- ,\Bi\) = ,h), the ratio of m's having \B\\ = b n is equal to 

b n /{bi + • • • + h) = bn/m- It is known that the number of w G NC(m) meeting conditions of 

• w has I blocks with sizes in a nonascending order of (&i, ■ ■ • , 6/), and 

• K{w) has m — I + 1 blocks with sizes in a nonascending order of (ci, ■ ■ ■ , c m _z + i), 

is m(m — /)!(/ — l)!//(6i, ■ • • , bi)/f(c\, ■ ■ ■ ,c m _/ + i) [43,10]. Thus, it is straightforward to see that (13) 
holds. 

Appendix IV 
Proof of Theorem 3 

A. Proof for jl m 

For m > 1 and 1 < i < m, let NC^\m) denote the set of noncrossing partitions in NC(m) such that 
the block containing 1 contains % as its largest element. Because of the noncrossing condition, a partition 
w £ NC^\m) can be decomposed into w = wiUzu2, where w\ G NC^\i), m<i G iVC({i+l, • • • , m}), 
and \w\ = \wi\ + |-ct72 1- It is clear that NC^\i) is in bijection with A^C(i — 1), so we can let w\ G 
NC({1,--- 

Consider a noncrossing partition w G NC^\m), and w = w\ U vj<i is the same as those described 
above. In its K-graph, there are \vo\\ + \ w2\ vertices, where vertices 1 to \wi\ correspond to blocks of 
vj\, and the (\w\\ + l)-th to the i\w\\ + |tU2|)-th vertices correspond to blocks of wi. Similarly, among 
the m — \w\\ — \w2\ + 1 cycles, cycles 1 to (i— \wi\) correspond to blocks of K(w\), and the remaining 
cycles correspond to blocks of K(w2)- 

To prove the result, the following observations are important. Among the cycles where vertex 1 is 
located, there is exactly one (the (i — \w\\ + l)-th cycle) corresponding to K(w2)\ all other cycles 
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correspond to K(w±). For any other vertex (other than vertex 1), either all its located cycles correspond 
to K(w\) or all the cycles correspond to K(w-i), but never some to K(w{) and some to K{w?). 

Suppose that w\ = {-Di, ■ ■ ■ , and 072 = {-D| roi |+i, • • • , Z)| ro |}. Owing to the above observa- 

tions, we can write (45) as 

M 

K -i N - m e n n v^ Sk ,b k ) 

6|wi|+i," - ,6|w| fc=|t»7i|+l s k £{B K (zv 2 ) 

Ci-|wi|+l," - >Cm-|w|+l 

|n7i| 

x E II V(c Sl ,bi) 

6i,->6|wil ' =1 s;e{BK(„ l) (') ; ' e - D i} 

= i^-^-™ e nq-Ki+i.^x En n v ^ h >) 

bi,Ci_| OT1 | + i 6 2 ,--- ,6|wii 1=1 si£{B K(zui) (t):t£Di} 

Ci ,Ci— 1 Wl 1 

M 

e n n ^(^,6*)- (53) 

+ - fe=|roi| + l s fc e{BK( ro2 )(<):iG-Dfc} 

c i-|i»ll+2>'" ,Cm-| w | + i 

Since NC(m) = U^iVC^m) and this is a disjoint union, we consider all 1 < i < m along with all 
w\ and W2, which expands (53) to 

61) (54) 

6i,c»_| TO1 |+i 

m |tJ7i 

xE^-^- (t - i)+i e e n n nc,® 

»=1 07ieJVC(i-l) 62,-,6|w!| «=1 si6{Sjr( wl )(t):teD,} 

roi={Di,- ,D| W1 |} Ci,"- ,Ci_| roi 

xjf-ijv-(™-o + i ^ e n n ^,6*)- 

vj 2 £NC(m—i) 6| wi |+i,— ,b| w | fc=| ro i|+l s fc€{£>K( ro2 )(t):t6-Dk} 

072 = {i3|wi|+l,— i-Dlwl} Ci-|wi|+2,— ,C m _|„| + i 

Note that there is a slight abuse of notational usage at the first line of (54): the summation variable 
Ci-|o7i|+i is indexed by variables % and tui that appear at the second line. However, as Cj_| roi |+i is a 
simply dummy variable, this abuse does not affect the result. 
Define 

M 

U m (c,b) := N- m+1 e e n n V ^'^ 



vueNC(m) 62,-,6|w| 1=1 Sie{Bft:( ro )(t):te-Bi} 

IC7 = {Bi,--- C2,-" ,Cm-|w|+l 



b±=b,c-i=c 
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Following the same steps in the proof of Theorem 1 for fi m , we have 

A m = P Km K^N- 1 V U m (c, b). 

K,N-^oo 
K/N^/3 



(55) 



b,c 



For each iV (with K/N = /?), let p!£P be a function given by 



P-nPfay) = U m (i,j), 



< X < — , < V < — , 

N - N K ~ K 



(56) 



and the limit of the sequence {p,^P(x, y)}^ =1 is jl m (x,y). Then, using the Riemann sum expression, 
(55) can be written as 

Am = P I \ fi m (x,y)dxdy. 
Jo Jo 

Let us represent the expression in (54) as Q(c,b) when the summation at the first line is discarded and 

b\ := b and Ci-\ Wl \ + i := c. We have 



P rAF 1 y]^(c,6)=Am- 

K /V — »rvi ' ■* 



K,N^oo 



(57) 



b,c 



Compare (55) and (57) and replace U m [c, b) as p,^P((c — 1)/N, (b — 1)/K), we can see 

n(c,b) = K-i N -^ (^,^ 

On the other hand, from the expression in (54), $7(c, b) can be also written as 

'c- 1 J-l' 



(58) 



ArV^) 



N ' K 



m—1 



i=i 



Cl 



^Iwil + l 



A 7 " 



K 



(59) 



where the sum YliLi m (54) is decomposed asl<i<m — 1 and i = m because tx7 2 at the third line 
of (54) is empty when % = m. Equating (58) and (59) and taking limits on both sides, we obtain the 
recursion in (16). 



B. Proof for f) m , k 

We prove this theorem by following the line of the proof for fi m . First, we derive the formula of f\ m ^. 
We multiply (54) by K and set b\ := k. Define 



M 

E m , k (c) := N-™+' Yl E II II V ^bi) 

vj£NC(m) & 2 ,-,b iH 1=1 s,e{B K (*,)(t):teB,} 
w={B u - ,B| w |}c 2> - ,c m -|„|+i 



bi=k,ci=c 



(60) 
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We have 



»7m,fc = I™ iV 1 V-E m ,fc( 



C ■ 



(61) 



For each AT, let be a function given by 



i — 1 i 

"aT^ty' 



and the limit of the sequence {j?^(x)}^ =1 is rj m ^{x). Then, (61) can be written as 



Vm,k 



JO 



Following the reasoning to obtaining (58) and (59), we have the equality of 



m— 1 



-(AT) / gi - 1 



N 



1 6i 



C7l|+1 



TV ' if 



+7V • A^ 1 V« W 



ci-1 



(62) 



By taking limits on both sides of (62), we acquire the recursion in (18). 



C. Proof for S mjk 

Suppose that we have S m ^(x) and fi m (x) such that 

S m ,k = 5 mt k(x)dx and fi m = (3 1 fl m (x)dx, 
Jo Jo 

where fi m {x) can be obtained from fi m (x, y) by integrating over y. Comparing formulas of fi m and b m ^ 
given in Theorem 1, we can see 

<W(x) = 0v k (x)fi m (x). (63) 
Integrating both sides of (63) over x G [0, 1), we obtain the relation of 5 m ^ and jl m (x,y) given in (17). 



Appendix V 
Proof of Theorem 4 

Here only the proof for /x m is given. The proofs for r\ m & and 5 m ^ can be proceeded in very similar 
ways. We start the proof from (45), and jJL m can be written as the sum of 



lim K- l N- m V Y E{ 

6l,'")&|w| Cl,— ,C m _| w |+l 
&l^-^6|w| Cl^-^C m _|„| + i 



K/N^f3 



n n 

k i=lsie{BK ( „)(t):teBi} 



(64) 
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over all noncrossing partitions w = {Bi,--- , B^i} £ NC(m), where the expectation is over the 
distribution of V(i,j). Since, for each j, {V(i,j) : 1 < i < N} is a realization of random process 
Zj(u,t), and Zj(u,t)'s are i.i.d. random processes for distinct j, the expectation in (64) is equal to 
M _ _ 

!'•, } I. (65) 



i=i 



H Z(u,c Si )\ =f[Mom z ( [j \ 

Jie{B K ^)(t):teBi} ) i=l Vs i e{Bif( ro) (t):t6B i } 



where the moments of Z(u,t) are employed to represent those of any Z^(u,t). Plugging (65) back to 
(64) and expressing the result as a Riemann sum, we obtain 



lim KMjV-M V" iv _ m+ | ro |_ 1 T-r < 

y jy , p Ci,-" ,c m _| w |+i i—1 



^i^MomJ [j {c Si } 

fe» Vs i e{B ft r (ro) (t):t6B i } J 



= p h f ■■■ [ n M ° m 4 u {* s ,}Wi---dx m _ H+1 , 

where the last equation is equal to the right-hand-side of (21) when summed over all noncrossing partitions 

w G NC(m). 

Appendix VI 
Proofs of Theorems 5 and 6 

A. Proof of Theorem 5 

The proofs for (24), (25) and (26) are very similar. Here only (24) is proved. We prove it by means 
of (21). We have 

Mom zlx (X^-i) = e| z(u,X Si )x\. (66) 

We use M(u), T(u) and W(u) to denote the real numbers that correspond to a sample point u G U for 
random variables M, T, and W, respectively. It is seen that, for any u eU, Y[ Sl e{B K (^)(t)-teB t } z ( u i X Si ) 
is nonzero (equal to M(u)l Bs l) if and only if either of the following conditions holds 

1) W{u) = 1 and maxjXi : I € \J teBi B K{w) (t)} < T{u) < 1, 

2) W(u) =0 and < T(u) < mm{x l : I G U teSi Bjf( w )W}. 

By considering distributions of T and W, events 1) and 2) above have probabilities 

1 , , 1 

- (1 - max^ ;i ) and - mm4 ;i , 



M 



respectively. Thus, the conditional moment in (66) can be reduced to 

E{Ml B 'l} ■ - (1 - vaas.X w . i + mmX m 
Plugging (67) back onto (21), we obtain (24). 



(67) 
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B. Proof of Theorem 6 

Only the formula of fi m is proved. Similarly to the proof of Theorem 5, ji m can be evaluated by 
means of (21). We can see W Sl e{B K (^){t)-teB t } z ( n ' -^0 * s nonzero (equal to M(u)^ B ^) if and only if 
z(u,X s .) / for all S{. This condition has the probability of P(u)\ Bi \. Thus, the conditional moment 
Mom z | X (A' ro; i) in (21) is reduced to E{M\ B ^}E{p\ B -\}. Plugging it back onto (21), we obtain 

M 

ft»= £ /? M n E { MlB<l } E { p|s<1 }' 

roeiVC(m) i=l 
n7={Bi,-,B M } 

which is equal to (27) by the fact: the number of Z-block noncrossing partitions having sizes of blocks 
in a nonascending order of b±, 62, ■ ■ ■ , h is equal to m(m — 1) ■ ■ ■ (m — I + 2)//(&i, 62, ■ ■ ■ , &j) [35]. 
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